Dirac and Normal Fermions in Graphite and Graphene: 
Implications to the Quantum Hall Effect 
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Spectral analysis of the Shubnikov de Haas (SdH) magnetoresistance oscillations and of the Quan- 
tum Hall Effect (QHE) measured in quasi-2D highly oriented pyrolytic graphite (HOPG) [Phys. Rev. 
Lett. 90, 156402 (2003)] reveals two types of carriers: normal (massive) electrons with Berry phase 
and Dirac-like (massless) holes with Berry phase tt. We demonstrate that recently reported integer- 
and semi-integer QHE for bi-layer and single-layer graphenes take place simultaneously in HOPG 
samples. 
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PACS numbers: 81.05.Uw, 71.20.-b 

Observations of the quantum Hall effect (QHE) and 
magnetic-field-driven metal-insulator transition in quasi- 
2D highly oriented pyrolytic graphite (HOPG) indi- 
cate that much in the physics of graphite had been missed 
in the past and had triggered a considerable interest in 
graphite once again. Analyzing the quantum de Haas 
van Alphen (dHvA) and Shubnikov de Haas (SdH) oscil- 
lations in bulk HOPG samples we experimentally proved 
that, besides the conventional electronic charge car- 
riers with the massive spectrum E — /2m, massless 
(2-|-l)-dimensional fermions (holes) with Dirac-like lin- 
ear spectrum E = ±v\p\ and non-trivial Berry phase tt 
do exist in graphite. It is likely that these Dirac-type 
carriers are responsible for the strongly correlated phe- 
nomena predicted by theory 0, d, . The results of Q 
were confirmed by the recent angle resolved photoemis- 
sion spectroscopy (ARPES) experiments performed on 
HOPG 0,1^, that provide an unambiguous experimental 
evidence for coexistence of Dirac-like holes and normal 
(massive) electrons located at H and K poins of the Bril- 
louin zone, respectively. 

Very recently, remarkable progress in the technology of 
ultrathin films consisting of one 

two [m or sev- 
eral graphite monolayers (graphenes) was achieved. 
Analysis of unconventional half-integer QHE and Berry 
phase TT of SdH oscillations in graphene led to the con- 
clusion that charge carriers in graphene have the nature 
of Dirac-like massless fermions [9|, [lOj. At the same time, 
bilayer graphene demonstrates a new type of the inte- 
ger QHE in which the last (zero-level) plateau is missing, 
and therefore the charge carriers are described by the 
chiral fermions having the Berry phase 27r and the en- 
ergy spectrum is represented by two touching parabolas 
E{p) = ±p72to [ill. 

In the present work we apply our filtering and phase 
analysis procedures Q to analyze the SdH oscillations 
and the QHE effect in HOPG aiming to clarify how 
these effects correlate with those studied in single and 



bilayer graphene systems. In particular, we demonstrate 
that normal and Dirac-like fermions seen in the bulk 
graphite are responsible respectively for the integer and 
semi- integer precursors of the QHE. To check the inde- 
pendence of these conclusions on the windowing and fre- 
quency filtration cut-off we used an alternative, insen- 
sitive to filtration procedure method of two-dimensional 
phase-frequency analysis, described in detail in 

We analyze the field dependencies of the basal-plane 
resistance RxxiB) and the Hall resistance Rxy{B) re- 
ported in Ref. [ij for the sample labeled as HOPG-3 
(noting that the HOPG-UC sample, similar to HOPG-3, 
has been studied in Ref. [§]). 

As seen from Fig. [U Rxy{B) demonstrates several 
QHE-like plateaus at high fields, also observed in vari- 
ous HOPG samples by other groups [H, 14 1. However, 
unlike the conventional QHE, Rxx{B) does not drop to 
zero at the positions of these plateaus. Only weak SdH 
oscillations /^Rxx{B) [3] are seen superimposed on al- 




FIG. 1: Field dependence of basal-plane magnetoresistance 
Rxx{B) and Hall resistance Rxy{B) measured for the HOPG- 
3 sample with magnetic field B applied along the sample c-axis 
[!]■ 
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most linear Rxx{B) dependence. This can be explained 
either by imperfection of the QHE itself or by the exis- 
tence of a background resistance of unknown origin. It 
is interesting to note that QHE with similar behavior of 
Rxx[B) has been reported for Bechgaard salts llSl lldl 
and for (Bii_2;Sb2:)2Te3 layered semiconductors [ITlllSt. 

Because quantum oscillations are periodic functions of 
the inverse field -B"^, the inverse- field spectral analysis 
is the appropriate tool to discriminate between different 
groups of carriers involved in oscillations. Fig. [^Ipresents 
the spectral intensity of SdH oscillations of the resis- 
tance Rxx{B~^) and the Hall conductance Gxy{B~^) = 
R~y{B^^), (in both cases polynomial background has 
been subtracted) together with analyzed in [3| spectral 
intensity of dHvA oscillations of magnetic susceptibility 
X(B-I), ix = dM/dB). 

Fig. [2] shows that two peaks in the susceptibility spec- 
trum observed at 4.68T and at 6.41T are also resolved in 
the spectrum of Gxy{B^^). The higher- frequency mag- 
netoresistance peak manifests itself as the right-shoulder 
structure of Rxx(B~^) spectrum. From comparative 
analysis of dHvA and SdH oscillations in Q we iden- 
tified these peaks as originating from normal (electrons) 
and Dirac (holes) carriers. Here, we re-affirm this con- 
clusion in a different way, namely by means of the com- 
parative analysis of longitudinal and Hall conductance 
oscillations, which is based on the separate study of the 
contributions from both types of carriers. 

Fig. [3] (curves a and b) shows the raw data of 
Gxy{B~^) and ARxx{B~^) (after background signal sub- 
traction) in which the normal electron contribution dom- 
inates in the magnetotransport quantum oscillations 
Extraction of the Dirac-like carriers signal from these 
data was the most challenging task that we performed 



by stop-band filtering of Gxy{B eliminating the fre- 
quency windows related to the electron spectral peak and 




FIG. 2: Spectral intensity of SdH oscillations of magnetoresis- 
tance \Rxx{Bo) \ and of Hall conductance \Gxy{Bo) \ in HOPG- 
3 sample in comparison with dHvA oscillations of susceptibil- 
ity |x(So)| in HOPG-UC sample % Peaks ei and h\ corre- 
spond to the normal electrons and Dirac holes, respectively 
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FIG. 3; Hall Conductance Gxy (solid line a) and l^a^x ~ 
—ARxx (solid line b) for HOPG-3 sample as a function of the 
normalized "filling factor" u — Bq/B (were Bq — 4. GST is the 
SdH oscillation frequency). Normal electrons give the dom- 
inant contribution in the quantum oscillations. The profile 
of Gxy{i^) is fitted by the integer QHE staircase; the minima 
of Aa^x correspond to the QHE plateaus. Hall conductance 
is normalized on Goxy ~ 28Q,~^ equal to the step between 
QHE plateaus. The results are compared with SdH oscil- 
lations in thick film of graphite (dashed line c) 13] (having 
approximately the same Bo = 4.65T) and with the QHE in 
bilayer graphene [Tl| obtained as a function of both magnetic 
field (dashed line d) {v = Bq/B, Bq — 27. ST) and concentra- 
tion (dotted line e) {u — n/no with no — 1.7 • W^^cm~^ at 



B = 20T) 
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FIG. 4: Solid lines a and b : the same as in Fig. [3] (with 
Bo = 6.41T, Goxy ~ 19f2~^) but after filtering elimination of 
the normal electrons contribution. The Gxy{i^) profile and 
positions of the minima in Aaxx ~ —ARxx are coherent 
with the Dirac QHE staircase. The results are compared 
with the Dirac QHE in graphene monolayer obtained both 
as a function of magnetic field {i/ — Bq/B with Bo = 6.67T') 
(dashed line c) [l3| and carrier concentration (v = n/no with 
no = 1.5 • lO^^cm"^ at B = 14r) (dotted line d) ^ 
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its harmonics (see Fig. [2]) and using the pass-band filter- 
ing of ARxx{B~^) with the corresponding hole frequency 
windowing. The resulting contributions from Dirac-like 
carriers are presented in Fig. 3] (curves a and b). 

Note that the sign of AR^x is reversed in order to 
present the oscillating part of the conductivity Aaxx, be- 
cause in our case of pxx ^ Pxy and Aa^x is related to 
ARxx as: 
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The situation is just opposite to the conventional QHE 
case with giant oscillations of Rxx{B) when pxx <C Pxy 
and a XX behaves in the same way as Rxx- Note also that 
the amplitude of SdH oscillations associated with Dirac- 
like fermions (Fig. 21 curve b) is three times smaller than 
that for normal carriers (Fig. [31 curve b). 

To facilitate the comparative analysis of the results 
obtained for various graphite (graphene) samples, we 
present all the data as a function of the normalized filling 
factor 1/ = Bq/B, where Bq corresponds to the SdH oscil- 
lation frequency. The conductance Gxy is normalized by 
Goxy, (— 28r2~^ for electrons and 19il^^ for holes) corre- 
sponding to the step between subsequent QHE plateaus. 

Fig. O (curve c) includes SdH oscillations reported for 
an 8pm thick HOPG sample taken from the inset of Fig. 
4 in [Tsll after the polynomial background signal subtrac- 
tion. The coincidence of the SdH oscillation frequency, 
phase and the asymmetric form of the signal with those 
measured in our HOPG sample suggests that the same 
Fermi surface electron pocket is responsible for SdH os- 
cillations in both samples. 

Before we proceed with the data analysis, we note that 
Landau Level (LL) quantization spectra for normal and 
Dirac-like carriers are different. In the normal carrier 
case the equidistant LLs En = {eh/m±c)B{n + 1/2) [l^ 
are separated by the gap Eq = ehB /2m±c from E = 
whereas in the Dirac-like case En — -^v yf2eKBTi[c [20| . 
and the Lowest Landau Level (LLL) is located exactly 
at Eq — 0. This leads to two important experimental 
consequences. 

(i) SdH oscillations of conductivity [1] 

Aa,,(B)~-A(B)cos[27r(:§-7-h5)] (2) 

(where A{B) is the non-oscillating amplitude) acquire 
the phase factor either 7 = 1/2 or 7 = for normal and 
Dirac carriers, respectively. The additional phase factor (5 
governed by the curvature of Fermi surface is quite small: 
|(5| < 1/8 (in 2D (5 = 0) and can be neglected. 

(ii) QHE plateaus occur either in integer (normal car- 
riers) or in semi-integer (Dirac fermions ) w ay and can be 
expressed via the phase factor 7 as 21 , S^l : 



where p = =f1 for electrons/holes, and gs is the (iso)-spin 
degeneracy factor. 

The phase factor 7 is defined for an arbitrary spec- 
trum E{p) by the quasiclassical quantization condition of 
the Fermi surface cross-section S{Ef) = {n + j)2TTeB / he 
(n » 1). This factor is uniquely related to the topologi- 
cal Berry phase $b = kir acquired by a fermion, moving 
around S{Ef) 7 equals either 1/2 for even k (nor- 
mal carriers) or for odd k (Dirac-like fermions). For 
example, the two-parabola spectrum of bilayer graphite 
E{p) = ±p^/2m_L is quantized in magnetic field as: 
En = ±{eh/m±^c)B^n{n — 1) This gives an un- 

conventional doubly degenerate LLL with Eq = Ei = Q 
but for higher n the system recovers the behavior of nor- 
mal carriers having 7 = 1/2 (with n — > n — 1 in p|)). 

We obtain the phase factors 7 and 5 by plotting the 
inverse field values at conductivity oscillation maxima 
(minima) as a function of their number (number - 1/2), 
as shown in Fig. [5l The linear extrapolation of data 
points to S"^ = unambiguously determines the phase 
factors: 7 = 1/2, 5 ~ —1/8 for data depicted in Fig. 
O (curve b) and 7 = 0, (5 = for data shown in Fig. [H 
(curve b). Thus, we identify here normal (3D) and Dirac- 
like (2D) charge carriers in agreement with our previous 
results obtained for HOPG-UC sample 0. 

Curves (a) in Figs. [H and [H provide evidence for the 
staircase behavior (quantization) of Gxy for both normal 
and Dirac-like fermions. As expected for QHE regime, 
the minima in longitudinal conductivity Uxx coincide 
with the plateau positions in Gxy 

In the case of normal carriers , Fig. [3l (curve a), the 
quantization of Gxy corresponds to the integer QHE, as 




Gx 
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(3) 
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FIG. 5: Values of for the normal (N) and Dirac-like 

(D)charge carriers at which nth Landau level crosses the 
Fermi level. These values are found as maxima (o) in the 
SdH oscillations of magnetoconductivity Aaxx{B). The min- 
ima positions (x) in Aaxx{B) are shifted by 1/2 to the left. 
Linear extrapolation of the data to — (solid lines) al- 
lows for the phase factor (— 7-1-5, Eq. (2)) definition: 7 = 1/2, 
(5 ~ — 1/8 for normal electrons and 7 = 0, 5 = for Dirac- fike 
holes. 
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follows from Eq. ([3]) for 7 = 1/2. In the same Fig. 
we presented the QHE staircases of axy, reported 
for bilayer graphite as function of both inverse field 
(Fig. [21 curve d) and carrier concentration n (Fig. [31 
curve e). The scales of and n were normalized to 
the corresponding oscillation periods of axx- The Hall 
conductivity was normalized to the step between sub- 
sequent plateaus in Uxy{B). The comparative analysis 
indicates that all three dependencies, i. e. Gxy[B^^) for 
HOPG, (Jxy{B~^) and axy{n) for bilayer graphene are 
equally well fit to the integer QHE theory with 7 = 1/2. 
We stress, however, that our data don't allow us to distin- 
guish between the conventional integer QHE with Berry 
phase and the chiral one with Berry phase 27r, proposed 



for bilayer graphite [Uj. This is because both models 



imply 7 = 1/2 and the same integer QHE staircase at 
LL n > 1. The only difference between them is the 
absence of zero level plateau for chiral fermions whose 
(non)existence we cannot verify at the moment. 

It is interesting to note that the absolute values of Hall 
steps A.p~y, estimated (similar to [ij) as (lOkfi/D)^^ for 
electrons and (ISkfi/D)"^ for holes are approximately 
equal to the two-spin degenerate QHE step 2e^//i ~ 
(12.9 kl7/n)~^ {gs = 2 in ([3])) and are twice as small 
as the two-spin two-valley degenerate QHE step with 
Os = 4, reported for graphene 0, and bilayer graphite 

To conclude, the results reported in this letter give new 
insight on the behavior of Dirac and normal fermions in 
graphite. Namely, we demonstrate the coexistence of 
QHE precursors associated with normal electrons and 
Dirac-like holes. To the best of our knowledge this is 
the first observation of the simultaneous occurrence of 
integer and semi-integer quantum Hall effects. 
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